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It is shown that the peak structure observed in angle-resolved photoelectron spectra of metallic 
deposits can only be unambiguously associated to single electronic states if the deposit has a two 
dimensional character (finite along one spatial direction). In one and zero dimensions the density of 
states shows peaks related to bunches of single electron states (the finer structure associated to the 
latter may not always be experimentally resolved). The characteristics of the peak structure strongly 
depend on the band dispersion in the energy region where they appear. Results for the density of 
states and photoemission yield for Ag crystallites on GaAs(llO) are presented and compared with 
experimental photoelectron spectra. 



The increasing availability of synchrotron sources is al- 
lowing the accumulation of experimental information on 
quantum size effects in metallic deposits on either metal 
or semiconductor surfaces-^ . This is for instance the 
case of silver deposits on Cu or GaAs substrates. Ag 
overlaycrs grow epitaxially on a Cu(lll) surface. In- 
stead, Ag deposition on the GaAs(llO) surface produce 
disordered films which, after annealing, give rise to Ag 
crystallites with a triangular cross section^. In both 
cases, Angularly Resolved Photoemission Electron Spec- 
troscopy (ARPES) show well defined peaks whose sepa- 
ration seems to be independent of energy, and decreases 
with the amount of deposit. 

In this work we show that these experimental data re- 
veal important dimensional effects which have not been 
discussed previously. Although a full understanding of 
the features of the photoelectron spectra would require 
a complete analysis of the specific electronic structure of 
a given deposit, some general trends can be derived and 
those are the object of the present work. We shall start 
by discussing the features of the free electron density of 
states in finite systems in two, one and zero dimensions. 
Then we shall present numerical results for the density of 
states and the photoemission yield calculated by means 
of a model suited to the case of silver deposits on the 
GaAs (110) surface. 

The density of states for free electrons in a. Lx M x N 
cluster is given by (atomic units will be used, h = 2m = 
e = l). 
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In the limit of an infinite solid the sums can be replaced 
by integrals giving the standard densities of states in one, 
two and three dimensions which show no oscillations^. 
Size effects can only be taken into account if these sums 
are not replaced by integrals. In two dimensions (a sys- 
tem finite in the I direction) m and n are good quantum 
numbers and the oscillations in the density of states are 



due to single electronic states. This interpretation was 
offered in^ to understand the origin of the quantum-well 
states observed in angle-resolved photoelectron spectra 
of Ag films epitaxially grown on Au(lll). This is not 
the case, however, in one and zero dimensions. In one 
dimension (a cluster finite in the I and m directions), the 
sums can be carried out in the limit of large energies by 
using Poisson formula^, leading to. 
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where En = E — {Trn/Ny and Jq is the Bessel function 
of the first kind and zero order. A constant density of 
states is recovered for an infinite solid. The accuracy of 
this expression is illustrated in Fig. 1, which shows the 
results obtained from Eq. (2) and those resulting from a 
numerical evaluation of Eq. (1) for the case of a 20 x 20 
cluster. It is noted that the peak structures are almost 
identical. However, whereas Eq. (2) gives an amplitude 
which decreases with energy, in the numerical results the 
amplitude increases. The reason is that, aiming to simu- 
late the almost constant broadening in the experimental 
spectra, the numerical calculations were carried out by 
adding a constant imaginary part to the energy, whereas 
the quadratic dispersion relation would have required an 
imaginary part dependent on the energy. 

The oscillations of Fig. 1 are not related to single 
electron states, but they are rather associated to several 
states each. It should be noted, however, that, depending 
upon experimental resolution, temperature and deposit 
dimensions, the finer structure associated to single elec- 
tron states may also be observed. The separation between 
these oscillations (peaks) decreases with the cluster size 
as 1/L^. On the other hand, peak separation, Ei — 
is not constant but rather paraboUic with increasing i. 
This is a consequence of the quadratic dispersion rela- 
tion characteristic of free electrons, and should therefore 
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strongljf depend on the band dispersion in the energy re- 
gion where the peaks He (see below). It can be checked 
that if one assumes a Unear dispersion relation (as can 
be found in some energy regions of actual band struc- 
tures of many metals^) the density of states in a L x L 
cluster shows peaks separated by a constant energy, and 
whose weight increases linearly with energy (in an infi- 
nite system such dispersion relation would give a density 
of states proportional to E). 
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FIG. 1. Numerical (continuous line) and approximated 
(Eq. (2)) results for the total density of electron states in 
a 20 X 20 cluster calculated within the free electron approx- 
imation. Results obtained with the approximated expression 
are plotted only for large values of the arguments of the Bessel 
function in Eq. (2). Both the density of states and the energy 
are given in atomic units {h = 2m = e = 1). 



As outlined above, the Ag 
GaAs(llO) have not a square 
rather have a triangular shape, 
that the density of states in a 
tions similar to those found in 
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crystallites formed on 
cross section but they 
It can be easily checked 
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a square. To illustrate 
case. The wave func- 
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the boundaries is given 
- sm{kxy) sin(fcyx)], C is 



a normalization constant. Note that the wave function 
within a square with infinite barriers at the boundaries 
can be expressed by either of the two terms in 'Q/. The 
dispersion relation is that of Eq. (1) with the restriction 
that I ^ m. The resulting density of states is essentially 
that of Fig. 1. 

These oscillations survive in the case of clusters in 
three dimensions. For instance, the density of states for 
a L X L X N cluster can be approximated by means of 
tecniques similar to those used in ID, giving. 
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As in the case of ID, this expression gives the correct 
result for an infinite system (only the first term in the 
rhs survives in this limit) and the agreement with the 
numerical results is again very satisfactory. The oscilla- 
tions are present but the structure is far more compli- 
cated. Let us consider the case N=2L, which roughly 
corresponds to the actual dimensions of the Ag crystal- 
lites considered below"'. The third term in the rhs of Eq. 
(3) approximately gives the same peak structure found 
in ID, whereas the second promotes a splitting of each 
peak into two. The picture which is roughly valid for N 
larger than L is the following: a density of states showing 
a gross structure similar to that found in ID with each 
peak splitted into approximately N/L peaks. This fine 
structure could not always be observed experimentally, 
depending upon resolution. As in ID, peak separation 
depends on energy. It is interesting to analyse the den- 
sity of states in a film (L = oo). In this case only the first 
and second terms survive. The first term accounts for 
the density of states of an infinite system in 3D, whereas 
the second gives oscillations of the type discussed above. 
This result suggests that angularly integrated photoe- 
mission spectra of metallic films would show oscillations 
similar to those already found in angle- resolved spectra-^ . 
Finally we note that a linear dispersion relation gives a 
peak structure with peaks separated by a constant energy 
and a weight proportional to . 

Before entering into a more detailed analysis of the 
experimental results for Ag deposits on GaAs(llO) sev- 
eral considerations arc in ordcr'^"^: i) The observed os- 
cillations sweap a rather large range of enegies, from the 
Fermi level down to around 4 eV; this is well below the 
top of the valence band of GaAs (less than 1 eV) . ii) The 
coupling of the crystallites with the substrate should be 
rather small due to the significant lattice mismatch be- 
tween the two solids" suggesting that these quantum well 
states weakly resonate with the valence states of GaAs. 
iii) The dispersion relation of Ag aroimd the Fermi energy 
is Hnear in k. iv) In the triangular cross section of the 
crytallites the Ag atoms are distributed on a rectangular 
lattice'' ((lie) plane of the fee lattice). 
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The above considerations, and the results for the free 
electron density of states in two and three dimensions, 
led us to carry out a calculation of the density states in 
triangular clusters of the square lattice by means of a 
tight binding model with one state per lattice site. This 
model can easily incorporate the rectangular lattice by 
taking different hopping integrals in the two directions of 
the square lattice; the three dimensional character of the 
crystallites can also be taken into account. Both, have 
been analysed and found not to change qualitatively the 
results herewith discussed, and their effects will be fully 
analysed elsewhere. Note that considering the fee lattice 
and the sp character of the conduction electrons in silver, 
would strongly limit the size of the clusters that could be 
reachable. The model, on the other hand, gives a disper- 
sion relation which, around the midband, is linear in k; 
note that in a metal (half filling) the midband coincides 
with the Fermi level. This overcomes the shortcoming 
of the free electron model discussed above which gives a 
quadratic dispersion relation for all energies. The trian- 
gular clusters considered here are bounded by the (01), 
(11) and (ll) fines, the former being the one assumed to 
be in contact with the substrate. This model cannot of 
course account for all the details of the geometry of the 
Ag crystafiites described in", but we think it takes into 
account their main features. Hereafter the unit of energy 
will be the hopping integral t (note that typical values 
of t adequate for normal metals would be in the range 
1.0-2.0 eV). 

The densities of states for triangular clusters having 51 
and 77 atoms in their respective bases arc shown in Fig. 
2 (the total number of atoms in a given triangular clus- 
ter with Nb atoms in its basis is {{Nb + l)/2)^). We first 
note that due to the linear character of the dispersion 
relation around the midband, the peaks are more regular 
in this energy region. Instead in the bottom and top of 
the band, where the dispersion relation is approximately 
quadratic, they are more irregularly distributed. Peak 
separation, on the other hand, is nearly 1.5 times larger 
in the smaller cluster. A fitting of the numerical results 
for clusters containing 576-9801 atoms gives a peak sepa- 
ration of Ei - Ei^i = l/(-0.24 + A^Nb) (with a correla- 
tion coefficient of 0.999) .It should be noted that although 
these features of the density of states qualitatively agree 
with the experimental data a detailed comparison would 
require a full calculation including s and d atomic or- 
bitals and clusters similar (in size and shape) to those 
grown in the experiment^'"'. 

This model has also been used to calculate the angle 
resolved photoemission intensity. The latter has been 
approximated by the following expression^. 



5{u -Eo-Wq + E,) (4) 



= 1,JV 



face of the photoemitted electron. In writing this expres- 
sion, the outgoing electron has been described by a plane 
wave exp{ikr) (with energy Eq = k"^ = {2tt I a)"^ + rn^) , 
where a is the interatomic distance in the square lat- 
tice). The electron states within the crystallites (en- 
ergy Ei) are linear combinations of atomic orbitals which 
have been approximated by (5-functions. A is the mean 
free path of the photoemitted electrons'^ and ds is the dis- 
tance of a given atom to the surface of the crystallite, or 
the distance to either the (11) or the (ll) lines (note that 
the cluster is assumed to be in touch with the substrate 
through the (01) line). 



0.5 



0.4 



0.3 



0.2 



CO 
LU 
I— 

CO 
O 

CO 

uu 

Q 



O.I 



o.o 



0.4 



0.3 



0.2 



0.1 



O.O 




-5.0 



-4.0 



-3.0 



-2.0 



-1 .0 



O.O 



where u is the photon energy, Wq the work function, and 
E and fc|| the energy and momentum parallel to the sur- 
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FIG. 2. Total density of states for two triangular clusters 
of the square lattice (defined by the (11), (11) and (01) lines) 
with 51 (a) and 77 (b) atoms in their bases, calculated within 
the tight binding approximation (one state per lattice site). 
Note that the density of states is symmetric with respect 
E = Q. The unit of energy is the hopping integral (see text). 

The results for the photoemission intensity perpendic- 
ular to the (11) or (11) lines (I = m) are shown in Fig. 
3. Figs. 3a and 3b shows the photoemitted intensity for 
LO — Wo = 18 and 12, respectively. The Fermi level is 
taken at the midband, A = 1.5a and (27r/a)^ = 10. It 
is clearly noted that although the relative weights of the 
peaks change, their positions remain unchanged. This is 
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a consequence of the reduced spatial dimensions of the 
cluster and is in agreement with the experimental data^'^. 
The effect of the emission direction is illustrated in Fig. 
3c, where the results for to — Wo = 12 and m = 0.75/ are 
plotted. We note important changes in the spectra, such 
as the dissappearance of some of the peaks, although the 
energies of the remaining peaks are not greatly modified. 
Again, these results qualitatively agree with the features 
observed in the experimental spectra. 
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FIG. 3. Angle resolved photoemission intensity for three 
sets of parameters: a) oj — Wo = 18 (in units of tlie flopping 
integral) and emission perpendicular to the (11) or (11) lines 
{I = m), b) same as a) for — Wo = 12, and c) same as b) for 
m = 0.75/. Energies are referred to Ep. The unit of energy 
is the hopping integral (see text). 

In summary, we have shown that the peak structure 
observed in photoelectron spectra of metaUic deposits 
can only be unambiguously associated to single electron 
states if the deposit has a two dimensional character 
(such as in angle-resolved spectra of films). In one and 
zero dimensions the density of states shows peaks related 
to many electronic states, showing a fine structure (re- 
lated to single electron states) which could eventually be 
resolved, depending upon experimental resolution, tem- 
perature and deposit dimensions. We have illustrated 



how the characteristics of the peak structure depend on 
the band dispersion in the energy region where they ap- 
pear. A tight binding calculation of the electron states in 
triangular clusters of the square lattice has been used to 
discuss some of the features of the photoelectron spectra 
obtained for Ag deposits on GaAs(llO). 
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